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Extraction of the Gluon Density, and Comments on "Hidden" Gluinos * 
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Calculations of deep inelastic processes (DIS) to next-to-next to leading order are discussed. Fitting ep experiment 
in the range 2.5 < < 230 GeV^, we find the coupling as{Mz) = 0.1163±0.0023. We also get the gluon density 
xG{x, ~ 10 GeV^) = 0.51 a;~'^ '**(l — x)^'^ , and negative evidence for the existence of light gluinos 



1. INTRODUCTION 

Deep inelastic scattering (DIS), in particular of 
electrons/muons on protons, constituted one of 
the first probes of hadron structure. The calcu- 
lation of QCD-induced scaling violations in the 
structure functions'^! yielded some of the earliest 
checks of the quark-gluon theory of hadron inter- 
actions, as well as providing the first two loop de- 
terminations of the strong coupling constant '^l. 

Let us set up some notation. Given the struc- 
ture function F2{x,Q^) in ep scattering, we split 
it into a singlet and a nonsinglet part, 

F2{x, Q2) = Fns{x, Q') + Fsix, g'). 

For the second we have to consider also the gluon 
structure function, Fc{x, Q^) = xG{x, Q^), G the 
gluon density because they mix. We project the 
moments, 



^^{n,Q^)= dxx"~'F,{n,Q'), i^NS,S,G 
Jq 

For NS the QCD NNLO evolution equation is 



-7;v'i(n)/2/3o 



1 + i?4^)a,(Q2)/47r + B^^\as{Q^) / Anf 



1 + Bi}'>as{Ql)/An + (a,(Q2)/4^)2 

The B may be written in tems of anomalous di- 
mensions, 7^^^(n) and Wilson coefficients, C^^\n) 
of order N: to NNLO we require iV = 0, 1, 2. (For 
the singlet equations, see ref. 3.) 

To compare the QCD predictions with exper- 
iment we need thus to evaluate 7'^^^(n), C^^^(n) 



(2), 



and invert the moments equations. This can be 
done with Altarelli-Parisi equations; but for this 
we would need the corresponding kernels, known 
only to NLO.W It is also possible to invert the 
equations analytically if the analytic form of the 
7'^^(n), C^^'(ri,) is known for all n (ref. 5 to LO, 
and ref. 2 to NLO). To NLO the calculations are 
complete[2^4.6]. nnLO we have the calculations 
of ref. 7 that provide us with the Wilson coeffi- 
cients; but the j^-^\n) are only known for a few 
values of n. Indeed, the calculations of Larin et 
al.I^l give those corresponding to NS scattering for 
n = 1; and the singlet and nonsinglet 7'^^-'(n) in 
electroproduction for n = 2, 4, 6, 8. Therefore, 
before comparing with experiment some work has 
to be done. For the nonsinglet case, see ref. 9; we 
next briefly describe the method followed by us in 
the singlet case, ref. 3 to where we send for more 
details. We also present here, for the first time, 
the ensuing determination of the gluon density, as 
well as a few comments on the (negative) impli- 
cations of our analysis for the existence of light 
gluinos. 

2. BERNSTEIN AVERAGES 

For a given value of only a limited number of 
experimental points, covering a partial range of 
values of x, are available, so one cannot simply 
use the moments equations. A method devised 
to deal with such a situation is that of averages 
with the (modified) Bernstein polynomials (modi- 
fied because, since only even moments are known, 
we have to consider polynomials in the variable 
x^); for details on the method, see refs. 3, 10, and 
work quoted there. We define these polynomials 
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for < n as 



Pnkix) 



2r(n+ I) 



„2fe 



r{k + ^)r{n-k + i) 



(l-x^) 



.2\n-k 



The Pnk{x) are positive and have a single max- 
imum located at Xnk ~ A;/n. They are concen- 
trated around this point, with a spread, Aa; 

nk 

1/n (accurate expressions for Xn^^ 

Axnk can be 

easily calculated, or looked for in ref. 3); and they 
are normalized to unity, Jq dxpnkix) = 1: so, 

the integral dx pnkix)(p{x) represents an av- 
erage of the function ip{x) in the region Xnk — 
^Axnk ^ X ^ Xnk + ^Axnk- The valucs of the 
function '•p{x) outside this interval contribute little 
to the integral, as p„k {x) decreases to zero quickly 
there. Finally, using the binomial expansion it fol- 
lows that the averages with the p„k of a function 
can be obtained in terms of its even moments: 



'nk = / dxp„kix)ip{x) 
■I 



2{n-k)\r{n+^) 



n—k 



i-iy 



r{k + i)r(n - + 1) § i\{n -k-iy. 



■ l^2k+2V 



and 1^21 = Jq dx x^'^ ip{x) . We will thus consider 
our experimental input to be given by averages 



2\ — 



dxp„fe(x)Fi^"P)(x,g2), (2.1) 



^^(exp) ^j^^ experimental structure function. I^^l The 
experimental points (and the theoretical fit) are 
shown in the figure. 
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NNLO WITH TUC 



X' = 79,2/(102-12) = 0.88 
a,(Mj) = 0,1 16J±0.0016 (ex 



-■0' 0- o ■f^--0.-.^{f_ 




3. NUMEROLOGY 

We present in Table 1 a compilation of the re- 
sults obtained with our calculations at LO, NLO 
and NNLO, with TMCs (target mass corrections) 
taken into account; the fit to the data itself is 
shown, for the NNLO calculation (with TMC) in 
the previous figure. Note the (small) jumps in 
the theoretical curves at the location of the mass 
threshold, =1^1; they occur because we have 
joined the theoretical formulas from n/ = 4 to 
ny = 5 at that point, using the method of ref. 12, 
which is not exact. 



Order 


A{nf = 4) 


a. (Ml) 


xVd.o.f. 


LO 


215 ± 73 


0.135 ± 0.007 


212/102 - 12 


NLO 


282 ± 40 


0.1175 ± 0.0027 


80.0/102 - 12 


NNLO 


283 ± 25 


0.1163 ±0.0016 


79.2/102 - 12 



Table 1 

We have 12 parameters: the twelve moments 



at the initial value, ^liin.Qf^) 



2, 4, 6, 8 and 



i = S,G, NS, minus one moment, ij,g{'2, Qq) which 
is deduced from /is (2, Qq) via the momentum sum 
rule; plus the QCD coupling parameter, A. The 
initial value is taken Qq = 8.75 GeV^, and we 
evolve to all other values of Q^, in the range 
2.5 GeV^ <Q^< 230 GeV^ with the QCD evo- 
lution equations; then construct the combinations 
entering the Bernstein averages, and fit the val- 
ues of these obtained from experiment. In Table 
1, only experimental {statistical) errors of the fit 
are shown; systematic (theoretical) errors will be 
discussed below. The NLO corrections are clearly 
seen in the fit: the x^/d.o.f. decreases from a lar- 
gish value of ^ 2.4 to a very good ~ 0.89. The fit 
is so good at this order that there is little room for 
improvement when going to NNLO; nevertheless, 
an improvement is seen. Not in the x^/d.o.f. ; but 
including NNLO corrections leads to a noticeable 
gain both in the quality of the determination of 
the coupling, and in the stability of the fits. 

Estimated systematic errors, originating from 
various sources, are shown for the NNLO calcu- 
lation in Table 2. No TMC means that we have 
neglect target mass corrections. The correspond- 
ing error is not included when evaluating the over- 
all theoretical error, since we do take into account 
TMCs in our central value. "Interp." refers to 
the theoretical errors inherent in the calculation of 
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the integrals in (2.1) giving the experimental aver- 
ages which arise because, for this calculation, it is 
necesary to interpolate the experimental points. 
We have used two difTcrcnt interpolation meth- 
ods, one assuming a hard Pomeron (refs. 13) but 
with independent fits for every Q^, which furnishes 
our central value, or using the MSRTgS,!!"*! that 
gives the estimated error. HT means that we have 
taken into account higher twists phenomenolog- 
ically, by adding, to /i7vs(^, Q^), the correction 
At^s(n, Q^) = n{aA^ /Q'^)iiNs{n, Q^). a is free pa- 
rameter whose fitted value is o = —0.202 ± 0.030, 
a very reasonable number. "Quark mass effect" 
means that we cut off the b quark threshold re- 
gion, instead of matching thruough quark thresh- 
olds; the variation in Table 2 takes into account 
also the variations due to the uncertainty in the 
mb mass. Qq to 12 GeV^ means that we take the 
input moments defined at this value of the momen- 
tum, i2i{n, Ql = 12 GeV^), i = S, G, NS instead 
of Qq = 8.75 GeV^ as was done to obtain the re- 
sults of Table 1. Finally, iVNNLO is the estimated 
effect of the (likely) larger sources of corrections 
of higher order in Ug . 



Source of error 


A(nf = 4, 3 loop) 


AA 


Aa,(M|) 


No TMC 


292 


9 


0.0006 


Interp. 


273 


10 


0.0007 


HT 


292 


9 


0.0006 


Quark mass effect 


299 


16 


0.0010 


to 12 GeV^ 


294 


11 


0.0007 


JVNNLO 


289 


6 


0.0004 



Table 2 

Composing quadratically all errors we find 

A{nf = 4, 3 loop) =283 ± 25 (stat.) ± 24 (syst.) 
=283 ± 35 MeV; 
af^°°p\Mz) =0.1163 ±0.0023. 



In Table 3 we compare our results to previ- 
ous determinations I^^' for as(Mf ), to the NNLO 
level (but excluding e+e" annihilations.) DIS 
means deep inelastic scattering, Bj stands for the 
Bjorken, GLS for the Gross-Llewellyn Smith sum 
rules. The xF^ result is that of ref. 9. 

The previously existing average, also taking 
into account NLO calculations, was q;s(M|) = 



0.118 ± 0.006; when including both our rcsidt and 
that of ref. 9 the new average and error become 

as(M|) = 0.1165 ±0.0016. 

We discuss briefly how the calculation could 
improve. Adding the values of 7(71) for a few more 
values of n would allow us to extend the range and 
increase the precision of our evaluations: alrcdy 
two more moments would probably decrease the 
error in A in some 30%. This is a difficult task, 
but it appears that Moch and Vermascren are on 
the way to it (see the paper by S. O. Moch, these 
Proceedings; it is even possible that the anlytic 
expresKioii be found for ~: (/0)- 



Process 


Average or 
range (GeV) 


a.(M|) 


T decays 


1 - (1.777)^ 


0.119 


Z — » hadrons 


(91.2)2 


0.124 


DIS; ly, Bj 


2.5 


r. 1 99+0.006 
'-'•-'-^^-0.009 


DIS; u, GLS 


3 


0.115 ± 0.006 


DIS; V, XF3 


5 - 100 


0.117 ± 0.010 


DIS: v/fip-F-, 


2..") - 2:!{) 


O.lKi:! ± ().()()23 


Average DIS 


2.5 - 230 


U.iiO»_o.o020 



Table 3 

4. THE GLUON STRUCTURE FUNCTION 

A spin-off from our results is that we also get mo- 
ments of the gluon structure function, in particu- 
lar at our starting value of Q"^ = Qo- 



Qi = 8.75 Qg = 12 

^g(2) = 0.242041 ± 0.003 0.2390 ± 0.0015 

Mg(4) = 0.0020 ± 0.00023 0.0020 ± 0.0055 

Hg{6) = 0.00119 ± 0.00050 0.00142 ± 0.0013 

^0(8) = (0.99 ± 0.86)10"^ (1.00 ± 0.84)10"^ 

The results, fairly stable, do not fix the gluon 
density; to get it, extra assumptions have to be 
made. Following the hard Pomeron model we take 
the functional form 

Faix, Q2) = xG{x, Q2) = Agx-^{1 - xf , 

fixing A = 0.44 (ref. 13) to avoid spureous minima. 
Fiting the moments gives 



4 



f = 8.75 GeV^ 
Ag = 0.52 



f = 12 GeV^ 
Ac = 0.51 



f = 8.75&12 
= 0.51 



V - 



.17 



X 



: 8.08 
= 2.4 



X 



8.11 
= 9.0 



the last if making a global fit. The value of the 
X^/d.o.f. = 9.0/(8 — 2) in this last case is quite 
reasonable: note that Aq and v should depend on 
Q^, so we are, in the last fit, obtaining an average. 

It is interesting to compare the value obtained 
for Ag with the predicition of the hard Pomeron 
model. We write'^^' 

Fsix,Q^) {el)AsiQ^)x-\ A = 0.44, 
MQ') = (e^)a.(g2)-'^+(i+^)Bs, Bs = const., 

Ag/Bs = [d+{l + A) - Dii(l + A)]/Di2(l + A) ~ 
4.82. Thus {el)Bs ~ 2.1 x lO"^^ the ballpark 
of the values found in ref. 13. 

5. COMMENTS ON "HIDDEN" GLUINOS 

It is recurrently suggested that "hidden" light 
gluinos could exist. I^^' They would alter the evolu- 
tion of as : in this respect, the agreement of our de- 
terminations with those made at smaller (r decay) 
or higher energy {Z decay) provides strong evi- 
dence against their existence. Direct negative evi- 
dence is obtained as follows. Let /d.o.f.(M) be 
the chi square per degree of freedom, taking into 
account only experimental points with < M^. 
If, at a given M = Mq a channel for the pro- 
duction of particles not taken into account in the 
analysis opened, then X^/d.o.f.(M) should jump 
at M = Mo, and would continue deteriorating for 
larger M, as we are fitting with wrong theoretical 
formulas: neglect of NLO corrections deteriorates 
the x^/d.o.f. substantially, and gluinos contribute 
already to LO. A plot of X^/d.o.f.(M) at differ- 
ent M is given in the figure (the full dots; the 
straight line is the ideal value X^ /d.o.f.(M) = 1). 
The increase of the X^ /d.o.f.(M) near M = mi, is 
clearly seen, due to our approximate treatment of 
the threshold; but, since we are using correct the- 
oretical formulas (that is, with rt/ = 5) above that 
value, /d.oI.{M) decreases as we get far from 
TOft. This shows the effectiveness of the method. 



1 .6 



1 .4 



xVdo' (M) 



C.8 - 

0.6 -•' 

Q 4 

C.2 - 











50 75 1 QO 1 25 1 50 1 75 200 225 250 



No other step, or deterioration is seen in the 
whole range. Now, the presence of a particle with 
mass Mq produces effects in the Bernstein aver- 
ages below ~ Mq because they involve inte- 
grals with the variable x. We are sensitive to ener- 
gies up to -EMax. ^ ^/{^ - x)Qlji^ /x ~ 30 GeV, 
and slightly above. So it follows that we can ex- 
clude gluinos with masses below 16 ^ 20 GeV. 
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